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ABSTRACT 

We present mid-infrared photometry of three very young brown dwarfs located 
in the p Ophiuchi star-forming region - GY5, GY11 and GY310 - obtained 
with the Subaru 8-meter telescope. All three sources were detected at 8.6 and 
11.7/mi, confirming the presence of significant mid- infrared excess arising from 
optically thick dusty disks. The spectral energy distributions of both GY310 
and GY11 exhibit strong evidence of flared disks; flat disks can be ruled out 
for these two brown dwarfs. The data for GY5 show large scatter, and are 
marginally consistent with both flared and flat configurations. Inner holes a few 
substellar radii in size are indicated in all three cases (and especially in GY11), 
in agreement with magnetospheric accretion models. Finally, our 9.7/im flux for 
GY310 implies silicate emission from small grains on the disk surface (though 
the data do not completely preclude larger grains with no silicate feature). Our 
results demonstrate that disks around young substellar objects are analogous 
to those girdling classical T Tauri stars, and exhibit a similar range of disk 
geometries and dust properties. 

Subject headings: stars: low mass, brown dwarfs - stars: pre-main-sequence - 
circumstellar matter - planetary systems - stars: formation 



1. Introduction 

Now that a large number of brown dwarfs are known, some with inferred masses close to 
those of extrasolar giant planets, their origin has become a subject of widespread interest and 
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fundamental importance. One key recent finding is that substellar objects harbor accretion 
disks and undergo a T Tauri phase similar to that of low-mass stars, possibly suggesting a 
common formation mechanism. The evidence for disks around brown dwarfs includes near- 
and mid-infrared excess (e.g., Tamura et al. 1998; Natta & Testi 2001; Natta et al. 2002; 
Jayawardhana et al. 2003), spectroscopic signatures of accretion (e.g., Jayawardhana, Mo- 
hanty & Basri 2002, 2003; Muzerolle et al. 2003) and millimeter-wave emission (Klein et 
al. 2003). In some cases, these disks appear to last for up to 10 million years (Mohanty, 
Jayawardhana & Barrado y Navascues 2003), a timescale comparable to that of disks sur- 
rounding low-mass stars. Moreover, there are now indications of jets/outflows from objects 
close to or below the substellar limit (Fernandez and Comeron 2001; Barrado y Navascues, 
Mohanty & Jayawardhana 2004), further strengthening the T Tauri analogy. 

Given the much lower mass and intrinsic luminosity of brown dwarfs relative to T Tauri 
stars, the preceding discoveries raise intriguing questions about how the detailed properties 
of brown dwarf disks compare to those of T Tauri disks. Observations in the mid-infrared, 
where optically thick disk emission dominates over the substellar photosphere, can provide 
important constraints on the disk models. Here we report ground-based photometry in the 
8.6, 9.7 and 11.7/im filters, spanning the 10/im silicate feature, of three young brown dwarfs 
in the nearby p Ophiuchi region (age<l Myr, d~150 pc), and use these data to explore their 
disk characteristics in comparison with models. 



2. Observations and Data Reduction 

All our photometry was obtained with the Cooled Mid-Infrared Camera and Spectro- 
graph (COMICS; Okamoto et al. 2003) on the Subaru telescope, on Mauna Kea. In imag- 
ing mode, the instrument offers a 42"x32" field of view, with a plate scale of 0.13"/pixel 
(diffraction-limited seeing at ~10/iin is 0.31", or ~2pix); we observed under seeing conditions 
of l"or better. We acquired data through 3 filters: a narrow-band filter centered on 8.6/im 
(width R3 OApm), and two broad-band filters centered on 9.7 and 11.7/im (widths « 1/im). 
[GY92]5 (ISO-Oph 30) and [GY92]310 (ISO-Oph 164) were observed on 2003 June 18 UT; 
[GY92]11 (ISO-Oph 33) was observed on 2003 July 16 UT. The flux standards HR5013 and 
HD 152880 were also observed on the former night, and HD 152880 alone on the latter. Data 
were obtained in standard chop-and-nod mode for science targets and calibrators. Time 
constraints only allowed observations at 8.6 and 11.7/im for GY11, and not at 9.7/im; the 
other two p Oph sources were observed in all three filters. Typical integration times were 
~ 1200s per filter for the science targets. 

Aperture photometry was performed on the dark-subtracted, flat-fielded and coadded 
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final images, using the apphot package in IRAF^ . Our final errors are < 0.05 mag in cali- 
bration for HR5013 and HD 152880, and ~0.15 mag in aperture photometry for GY5 and 
GY310, resulting in a total error of ~15% in the inferred flux in all filters for the latter two 
objects. For the faintest source GY11, we conservatively estimate an uncertainty of ~0.25 
mag in our photometry, yielding a final flux error of ~25% at both 8.6 and 11.7/xm. 

GY5 is observed but not detected at 9.7/xm, in either the individual nod frames or in the 
coadded frames at a single nod position. Under the circumstances, we derive a 3a detection 
upper limit of 35 mJy from the individual frames. This value is consistent with our detected 
9.7/Um flux of ~35 mJy in GY310, observed for nearly the same length of time: the latter 
source is just visible at the 2-3<r level in the individual nod frames (but clearly detected at 
~6cr in the co-added frames). Our Subaru fluxes are listed in Table 1. In our analysis of the 
spectral energy distributions (SEDs), we have also included the ISO 6.7 and 14.3 /iin fluxes 
derived for these sources by Bontemps et al. (2001) with an angular resolution of 6 arcsec, 
and the ISO 3.6, 4.5 and 6.0 /im fluxes for GY5 and 11 from Comeron et al. (1998) with 
3-arcsec resolution; these all have errors > 10%. 



3. Stellar Parameters 

We are interested here in the luminosity, mass and radius of our targets. Mass is 
pertinent for establishing whether these are bona-fide brown dwarfs, while luminosity and 
(to a lesser extent) mass and radius are important for modelling their disk properties (§4). 
Our stellar parameters for GY5, 11 and 310 are adopted from Natta et al. (2002, hereafter 
NTC02); the reader is referred to that work for a detailed discussion of the methodology. 
The derived quantities are listed in Table 1. NTC02 estimate errors of ±100K in T e ff , 
20-30% in Lb i , and 20-30% in mass (the h^oi errors derive mainly from those in Ay , and 
are likely somewhat higher in the faintest, heavily reddened objects; see discussion of GY11 
in §5). They find a spectral type of M6, T e// of 2700K and mass of 60±20 Mj for both GY5 
and 310, and M8.5, 2400K and 10±2 Mj for GY11. The inferred bolometric luminosities 
for GY5 and 310 are very similar, 0.07 and 0.09 L respectively; GY11 appears much less 
luminous at ~0.008 L . These values put all three bodies below the brown dwarf boundary 
at ~80 Mj ; indeed, GY11 seems to lie close to the deuterium-burning limit of ~12 Mj (as 
also argued by Testi et al. 2002). 

For comparison, we note that Luhman & Rieke (1999; hereafter LR99) have also inves- 
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tigated the properties of p Oph members, using K-band spectra and an analysis technique 
differing from that of NTC02. For GY310, their inferred and T e ff , and hence mass from 
theoretical tracks, are almost exactly those derived by NTC02. The brown dwarf status 
of this object thus seems secure. For GY5, their T e jj and (2900K and 0.11 L ) are 
slightly higher than NTC02's values, implying a somewhat larger mass of ~ 80 Mj . We 
cannot be confident therefore that this is a bona-fide brown dwarf; nevertheless, this higher 
mass estimate still implies an object very close to the substellar boundary. Moreover, using 
LR99's estimates does not significantly alter our disk models for this object or affect our 
conclusions regarding its disk properties. 

The situation is somewhat more complicated for GY11. LR99 find this object to be 
simultaneously much hotter (T e// ~ 2800K) and fainter (L bol ~ 0.003 L ) than NTC02 
derive. Results similar to LR99's are obtained by Wilking, Greene & Meyer (1999; hereafter 
WGM99) as well; like LR99, WGM99 also employ K-band spectra. Theoretical evolutionary 
tracks (e.g., D'Antona & Mazzitelli 1997; Chabrier et al. 2000) then imply a substantially 
higher mass, 30-40 Mj , than NTC02's value. However, they also then imply an age of ~10 
Myr (see comparisons to evolutionary tracks in LR99 and WGM99). Now, the infrared excess 
and veiling of GY11, signifying large amounts of circum- (sub) stellar material, confirm it to 
be a true member of the young p Oph star-forming region (WGM99), which has an estimated 
age < 1 Myr. Thus, the LR99 (and WGM99) values for GY11 appear suspect. As WGM99 
point out, it is quite possible that veiling by dust causes its T e /j from a K-band analysis to be 
significantly overestimated. Moreover, the H 2 bands used by LR99 and WMG99 possibly 
saturate with later M type, which can also result in a spuriously high T e jj determination 
(LR99). Spectral type and T e ff errors will lead to erroneous inferred luminosity as well. 
The estimates by NTC02 are less susceptible to these particular uncertainties: their J- and 
H-b&nd analysis of the overall spectral shape depends less on H 2 band-strengths alone, and 
is plausibly less affected by excess emission from dusty disks than K-band studies. These 
considerations, and the fact that NTC02's luminosity and T e ff are in fact consistent with an 
age < 1 Myr, suggest that the NTC02 parameters we adopt for GY11 are more reasonable 
and appropriate. In any case, even the higher T e ff of LR99 implies a substellar mass for this 
object, whether one adopts their L bo ; or the NTC02 value (~ 40 Mj in both cases; also see 
§5). 



4. Disk Models 



To investigate the disk SEDs of our targets, we employ passive disk models that only 
consider the effect of reprocessed stellar irradiation and ignore viscous accretion heating; 
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this is justified for the very low accretion rates, ~ 10~ 9 -10~ 12 M yr -1 , inferred so far 
in the substellar regime (Muzerolle et al. 2003). Our models are based on the 2-layer 
description of Chiang and Goldreich (1997; hereafter CG97) as modified by Dullemond et 
al. (2001); scattering is included in an approximate fashion, following Dullemond & Natta 
(2003). The disk geometry is assumed to be either flared or flat. In the former case, 
consistent with hydrostatic equilibrium for well-mixed gas and dust (Kenyon & Hartmann 
1987), the vertical temperature structure at any radius is described by two components: the 
temperature in the interior of the disk and that at the disk surface, which is an optically 
thin layer directly exposed to stellar radiation. The flat, geometrically thin case is modelled 
simply as an optically thick blackbody disk which has no silicate emission. This situation 
may correspond to complete settling of large grains to the disk midplane (e.g., D'Alessio et 
al. 2001). We assume a distance of 150pc to the p Oph region, consistent with the latest 
estimate by de Zeeuw et al. (1999). 

The disks are scaled-down version of typical T Tauri disks; detailed parameters are given 
in NTC02 (see also Fig.l caption). The gas and dust are assumed to be well-mixed in the 
flared models. As discussed in NTC02, the emission in the near- and mid-IR depends only 
on a few disk parameters, namely the disk inclination and inner radius and the properties 
of grains on the disk surface, for which we adopt a mixture of graphite (containg 30% of 
the cosmic C abundance) and astronomical silicates (containg all cosmic Si abundance) with 
grain radius a. Two distinct particle distributions are examined: (i) a min = 0.2/im and a max 
= 3.6yum; and (it) a min =3fim and a max =12fim. In both cases, the size distribution goes as 
dn/da oc a~ 3 5 . The two situations allow us to explore the effect of grain growth on the 
SEDs. For the calculations presented here, the surface grains properties are not crucial as 
long as a remains smaller than a few microns. As for stellar parameters, the most important 
one is L bot ; for flared disks, there is also a weak dependence of the flaring angle on the stellar 
mass and radius (CG97). 



5. Results and Discussion 

SED model fits to our data for GY310, 11 and 5 are illustrated in Fig. 1. We plot both 
our new Subaru data at 8.6, 9.7 and 11.7/im as well as additional flux measurements at shorter 
and longer wavelengths, mainly from ISO. The flux-calibrated near-infrared (0.85-2.45/xm) 
stellar spectra from NTC02 are also shown. All the observed fluxes have been corrected for 
reddening following NTC02. For each object, we plot two flared disk models, one with a 
distribution of small grains and the other with larger grains; a flat disk (blackbody) model; 
and a model stellar spectrum (Allard et al. 2001; see NTC02) depicting the contribution 
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of the underlying stellar photosphere. The models incorporate inner holes of various sizes. 
The effect of an inner hole in a flared disk is twofold. On the one hand, it depresses the 
near-IR emission, by getting rid of the innermost hot dust; on the other hand, it increases the 
fraction of the stellar surface "seen" by the disk (from 1/2 to 1, in the extreme cases R{ — R* 
(no hole) and Ri 3> R* respectively), and therefore the emission at longer wavelengths. All 
three targets exhibit considerable mid-IR excess above the stellar continuum: while their 
earlier ISO detections left open the possibility of source confusion due to ISO's coarse spatial 
resolution, the Subaru detections confirm that these objects harbor circum- (sub) stellar disks. 
The salient results from our present work are as follows: (1) both GY310 and GY11 clearly 
possess flared disks - flat disks can be confidently ruled out in their case, (2) the 9.7/xm flux 
for GY310 indicates possible silicate emission from small surface grains (though larger grains 
with no silicate emission cannot be completely excluded), and (3) inner holes of size ~ 3-7 
stellar radii seem necessary to fit the data, analogous to the situation in classical T Tauri 
stars. We now briefly discuss each source. 

GY 310: Our best fit for this object requires a flared disk, tilted 60° from face-on and 
with an inner- hole radius of Ri ~ 7R*. Inner holes significantly smaller than this produce too 
much emission at 6.7/xm, while smaller inclinations overestimate the mid-IR flux. Notice that 
a flat disk model is completely inconsistent with the data longward of 6.7^m; the predicted 
fluxes are simply too low. Moreover, all the fluxes from 6.7 to 14.3 jum agree remarkably well 
with the flared disk SED that includes lOyum silicate emission from small surface grains. A 
silicate feature is routinely observed in disks around T Tauri stars (e.g., Cohen & Witteborn 
1985; Natta, Meyer & Beckwith 2000). If our detection is confirmed, it would be the first 
evidence of such emission from a brown dwarf disk. Nevertheless, we caution that our present 
data cannot unequivocally rule out a flared disk with larger grains and no silicate feature: 
our 9.7/um flux remains marginally consistent with the latter scenario. Measurements with 
smaller uncertainties, such as Spitzer low-resolution spectra, can settle this issue conclusively. 

GY 11: Our best-fit model for this brown dwarf is a flared disk seen face-on, with an 
inner hole of size Ri ~ 3R*. Much smaller inner holes produce unacceptably low mid-IR 
fluxes. Similarly, a flat disk produces emission wholly inadequate to match the steeply rising 
mid-IR SED. However, in the absence of 9.7 /im data, we are not sensitive to the presence of 
lOfim silicate emission, and thus cannot distinguish between flared disks with small surface 
grains and those with somewhat larger particles. We note that the Lbd we adopt for GY11 
(0.012 L ), to obtain an adequate SED fit, is moderately higher than derived by NTC02 
(0.008 L ). The NTC02 value produces mid-IR fluxes somewhat lower than observed (fit 
not shown); since we already maximize the mid-IR emission by adopting a face-on disk 
and invoking an inner hole~ 3i?*, higher luminosity seems the only remaining option in 
the context of our models. The adopted here is admissible within NTC02's Ay errors 
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for this heavily reddened object (our L bot amounts to A v ?^8.5 instead of NTC02's 7.5±1). 
The T e ff from synthetic spectral fits assuming our new Ay is 2500K, very close to NTC02's 
2400K; the new L bo i and T e ff also remain consistent with an age < 1 Myr. The mass now 
inferred from the evolutionary tracks is ~ 20 Mj ; while higher than NTC02's estimate of 
10 Mj , it still implies that GY11 is a very low- mass brown dwarf. 

GY 5: The large scatter in the GY 5 data make any fit provisional at best. However, 
if we disregard the ISO 3.6/zm flux, which is 50% larger than the ground-based value in the 
comparable L' filter (from Comeron et al. 1998 and Jayawardhana et al. 2003) and assume 
that the true fluxes at 6.7 and 14.3^m are at the lower end permissible by the ISO error 
bars, then a flat disk inclined at ~ 65° to face-on, with an inner hole Ri ~ 5/?*, provides a 
plausible fit (especially since our 9.7/zm Subaru estimate is only a flux upper limit). However, 
we cannot discount a flared disk seen at large inclination. 

The disk parameters presented here for the three p Oph sources are not unique: there 
are plenty of trade-offs between factors such as disk inclination, inner-hole size, grain albedos 
and sizes, and degree of dust settling. For example, our models adopt only the two extremes 
of grain settling: well-mixed gas and dust leading to full disk flaring and complete settling 
resulting in flat disks. If there is some, but not complete, dust settling in GY310 for instance, 
its moderate mid-IR excess may be fit by a disk model with a lower inclination. 

Nonetheless, the SEDs we observe strongly suggest certain qualitative features that are 
relatively independent of the particular disk parameters adopted. First, the presence of 
flared disks in GY310 and 11 is indubitable: no flat disk, whatever its inner-hole radius or 
inclination, can reproduce the observed increase in their flux from 6.7 to 14.3 fim. Apai 
et al. (2002) carried out observations similar to ours for a single young object at the stel- 
lar /substellar edge, ChaHa2, and found that a flat disk with no silicate feature provides the 
best fit to its SED. Our results indicate that, while such disks might indeed be present in 
some cases (e.g., perhaps GY5), this does not signify a fundamental change in disk properties 
upon crossing from the stellar to the substellar regime: brown dwarfs can have flared disks 
too, as well as (tentatively) silicate emission, and simply present a range of disk geometries 
just as T Tauri stars do (e.g., Miyake & Nakagawa 1995). Second, inner holes of a few sub- 
stellar radii may be common in brown dwarf disks - our best fits to all three objects require 
such a configuration. Especially for GY11, it is hard to find an alternative; even assuming a 
face-on disk and a (sub) stellar luminosity higher than previously estimated, some mechanism 
is required to further raise the mid-IR fluxes to match the data, and a central hole seems the 
only viable possibility (see also Liu et al. 2003). Indeed, magnetospheric accretion models 
predict the formation of inner holes of just these sizes in both classical T Tauri stars and 
accreting brown dwarfs, due to disk truncation by interaction with the (sub) stellar magnetic 
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field (Shu et al. 1994; Muzerolle et al. 2003). In summary, therefore, our study suggests 
that brown dwarf disk properties are qualitatively similar to those of higher mass T Tauri 
stars. A more detailed understanding of grain characteristics, dust settling and geometry in 
brown dwarf disks will be possible with mid-IR Spitzer observations. 
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A (yu-m) 



Fig. 1. — Top panel: Observed fluxes and model predictions for GY310 (ISO#164). Blue 
dots from this paper, red dots from Bontemps et al. (2000); black curve shows the near-IR 
low-resolution spectrum from NTC02. All observations are corrected for reddening (Ay =6 
mag; NTC02). The magenta curve shows the prediction of a star + flared disk model with 
small grains; the blue curve a star + flared disk with large grains; the green curve is the 
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SED of a flat blackbody disk. The photospheric flux for the adopted substellar parameters 
is in yellow. The flared disks have Ri ~ 7i?*, inclination of 60°. The flat disk has Ri ~ 5/?*, 
same inclination. Middle Panel: Same for GY11 (ISO#033). The red squares are ISO fluxes 
from Comeron et al. (1998). The adopted stellar parameters (see §5) are =O.O12L , 
Ay =8.5 mag. All disks have Ri ~ 3i?* and are seen face-on. Bottom Panel: Same for GY5 
(ISO#030). The black diamond is an L' flux ground-based measurement from Comeron et 
al. (1998). The flared disks have Ri ~ inclination ~ 60°. The flat disk has Ri ~ 5R* 
and inclination ~ 65°. For all 3 objects, the other disk parameters are: disk outer radius of 
R D = lxlO 15 cm (67 AU); disk mass of M D = 0.03 M*; disk surface density profile S oc R~ l . 
The dust in the disk interior is assumed to have an opacity k = 0.01(A/1.3 mm) -1 cm 2 g _1 
(Beckwith et al. 1990). Overall results in this figure do not depend on the exact values of 
these parameters. 
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Table 1. Stellar Parameters and Mid- Infrared Fluxes 



name 


SpT a 


T a 

(K) 


Lfooi a 
(Lo) 


M a 
(Mj) 


A v a 
(mag) 


6.7/xm b 
(mJy) 


8.6//m c 
(mJy) 


9.7 (im° 
(mJy) 


11.7^m c 

(mJy) 


14.3^m b 
(mJy) 


GY5 


M6 


2700 


0.07 


60 


3.0 


28±7 


19±3 


<35 


24±4 


33±14 


GY11 


M8.5 


2400 


0.008 


10 


7.0 


9±3 


11±3 




21±5 


17±3 






2500 


0.012 


20 


8.5 












GY310 


M6 


2700 


0.09 


60 


6.0 


20±2 


24±4 


35±5 


39±6 


39±6 



a Spectral type, T e ff , Lboi , mass and Ay from NTC02. For GY11, we also provide, in the second 
row, the moderately different parameters suggested by our SED fits. 

b 6.7 and 14.3^m ISO fluxes from Bontemps et al. (2001). 
c 8.6, 9.7 and 11.7 fim fluxes observed by us at Subaru. 



